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Abstract 

Using improved mean field and strong coupling expansions we re-analyse the bulk phase 
diagram of the fundamental-adjoint action of the SU(2) Lattice Gauge Theory. We find that 
the qualitative features of the bulk phase diagram are robust and unchanged by the inclusion of 
higher order terms. On the other hand, some of the quantitative features, such as the location 
of the endpoint of the line of bulk phase transitions, seem to be strongly dependent on the 
higher order terms of the strong coupling expansion. 
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Lattice regularization of a continuum action is not unique. For non-abelian gauge theories 
the Wilson action jl]] is the most popular one, but other actions have been studied in the 
literature. In particular, lately there has been a resurgence of interest in the Bhanot-Creutz 
action |Q] 

S = E{/ ? /( 1 " \ttfU p ) + (3 a (l - \Tr A U p )}. (1) 

p 

Here F and A denote the fundamental and adjoint representations respectively. The Wilson 
action is a special case of (|]), corresponding to (3 a = 0. The action (|l]) was first studied by 
Bhanot and Creutz for SU(2) gauge theories in order to understand the bulk phase transition 
found in numerical studies of some non-abelian gauge theories (SO(3), SU(4) etc.) and the role 
they play in the physics of confinement. They found a line of first order transition in the (3f - 
[3 a plane (see Fig.l) that ended at a finite (3 a . Since the location of the peak in the plaquette 
susceptibility for the Wilson action corresponds to the interception of the extrapolation of this 
line with the (5 a = axis, it has been considered as a possible source for the observed crossover in 
the string tension. It is thus a possibly important key in our understanding of the confinement 
phenomenon. 

However, recent ]|, ||| finite temperature investigations of this action have cast some doubt 
on the nature of the phase transition line seen in Q . It was found that switching on a nonzero 
f3 a , the known finite temperature phase transition of the Wilson action becomes a line and joins 
the above mentioned bulk transition line. Moreover, the order of the deconfinement transition 
changes from second to first order at (3 a > 1.5. No indications of two separate transitions were 
found at any j3 a . After considering various possibilities, it was concluded Qthat the transition 
line is not a bulk one, but the deconfinement transition line. Since the study in || was done on 
small lattices, which were at relatively high temperature, such a misinterpretation is possible. 

A finite size scaling-based analysis of the bulk transition has also been done @ to determine 
the nature of the phase transition. It was found that the line of the phase transition ends at a 
somewhat higher (3 a than what was found in Q. Simulations at a still higher value of (3 a (=1.5) 
suggested the presence of a 1st order bulk phase transition, but the Polyakov line, which is the 
order parameter of the finite temperature phase transition, was also found to jump across this 
transition at j3 a = 1.5, thus further making a distinction between the zero temperature bulk 
transition and the finite temperature deconfinement transition very difficult. 

Numerical simulations thus seem to give conflicting signals. The deconfinement order 
parameter acquires non-zero vacuum expectation value at these transitions for all (3 a . However, 
the shift of the transition point with N% - the temporal lattice size - decreases and becomes 
negligibly small. The latter is suggestive of a bulk transition unless this behaviour changes on 
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very large lattices. In the absence of such large lattice simulations it may be instructive to look 
for guidance by conventional analytical techniques. 

The existence of a deconfinement phase transition for (5 a = has been rigorously proven 
in the strong coupling limit. In || a leading order strong coupling analysis of the action 
in (|l]) was done at finite temperature. It yielded a deconfinement phase transition line in the 
(/3f,P a ) plane, along with a change in order of the phase transition for large p a , as seen in ||. 

A mean-field analysis of ([j]) at zero temperature was done in Ref || and the results of (2| 
were supported. In view of the results of ||, however, it seems to be necessary to reexamine 
the results of Ref In this note, we have attempted to check how stable the results of Ref 
H are, both qualitatively and quantitatively, by improving and extending their work to higher 
orders in a consistent manner. We find that our results still predict qualitatively the same bulk 
phase diagram as in (|| but the location of the endpoint of the transition line is very sensitive 
to higher orders and cannot be precisely obtained by our study. 

The mean field analysis proceeds by writing the partition function 

Z(/3 f ,p a )= fUdUteM-S) (2) 
J i 

in the axial gauge by fixing all the links in time direction equal to 1. The SU(2) elements are 
parametrized as U = u° + iu.a, for real numbers (uq, u) satisfying Uq + u\ + u 2 . + = 1. The 



standard technique M, 10 1 of Fourier transforming the measure and the action gives 



Z((3 f ,(3 a ) = Jj[dv l j "l[J^eMS(/3f,M + - afvf - a,.v,)) (3) 

where 

w(a) = ln J cttlexp(aV + a.u). (4) 
One then looks for translationally invariant saddle points of (||) of the form ||] 

vi = (u,0),a, = (a,0). (5) 

For (3f < 2 one gets the stable solution v = 0, a = [Q. An expansion of the free 
energy around this saddle point gives just the strong coupling expansion. The strong coupling 
expansion for the Bhanot-Creutz action has been obtained upto terms of order sixteen by 
Dashen et al [O]. In our notation the series for free energy as obtained by them is 

Fi{p f , P a ) = 61n6 + 4(46? + 96^) + 36(4^° + 9^°) 

+36(12^% + 2Ab\b\ + 276^ - 224&} 2 - 13686?^ - 13596^ 2 ) 
+4(4326f6|6 3 + 405^% + 166| + 25fe|) (6) 
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and the characters 60 and bj, j=l,2,3,4, are evaluated as 
dO . , 9 r „ 9,2 



bo 



PITT fjfj ff / fj O s 

j Q -^sin 2 %xp(-/3 / (l-cos|)--/3 a (l-cos0)), (7) 
h i = JjTW l ^sin^sm^±^exp(-^(l-cos^)-^ a (l-cos0)). (8) 

The characters are now expanded in a series of /3f and (5 a to rewrite Fj explicitly as an expansion 
in powers of (if, j3 a . We checked that for our region of interest, the difference between eqn(|6|) 
and the explicit series is negligible. Also the convergence of the strong coupling series was 
checked by looking at the series terminated at differet orders. It was found that the strong 
coupling series converges very slowly near the transition line, specially for smaller (3 a values. 

For regions of large (5 a and (5f there are stable solutions of (0) for nonzero v,a satisfying 
the equations || 

» = 'M < 9 > 

a = /3 f 2v(l + (d-2)v 2 )+p a ^-v 3 (l + {d-2)v A ), (10) 

where I2, I\ are the modified Bessel functions; d, the space-time dimensionality, will be 4 in all 
our calculations. 

Expanding near this saddle point, the partition function, or equivalently, the free energy 
per site, is 

Fii(J3f,P a ) = In Z(J3 f ,p a ) = F tree + F^ loop + F 2 _ loop + ... (11) 
The first two terms of this series were calculated in Ref S and are 

F tree = 3f3 f (v 2 + v i -2)+4 f 3 a (v i + v 8 - 2) + 3(ln^M - <«,), (12) 

Fi-ioop = ^(ln2 + 31n|--ln(/3 / + ^ 2 )) -31n(/3/(l + 3v 2 ) + ^(3 a v 2 (I + 3v 4 

+3K X (13) 



where 



with 



and 



f w d 4 p 

Kl = ~J (- 2 7t) 4 ~ Cl cosp ° ~ ^2 (cos pi + cos P2 + cosp 3 )) (14) 

C x = (l + 3c)- 1 ,C 2 = (3 + -)- 1 (15) 

c 

_ 3(3 f v 2 + 8(3 a v 6 
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A comparison of magnitude of Ft ree and F\-i oap near the phase transition line shows that 
they are of the same order of magnitude, and so care should be taken in using (pd|). The loop 
expansion in ( |TTj ) is justified by saying that the higher loop corrections are suppressed by inverse 
powers of a J10|J . However, the lower loop terms might have a large contribution from higher 
powers of 1/a. Reexpanding Ft ree + Fi_i oop in ( pd| ) we found that the term proportional to 1/a 
comes with a coefficient ~ O(10). So for consistency, one should explicitly write the series as 
a series in inverse powers of a. Both due to the sensitivity of the location of the endpoint of 
the phase transition line, and to check the convergence of the series in 1/a, it is necessary to 



evaluate the series up to the term proportional to 1/a. This can be done [14] by a comparison 
with the weak coupling perturbation series. Using the weak coupling series of Ref [11], one gets 
the following expression for free energy around this saddle point : 



2 4- — q r I 

F H = -^fa lna + c + ^ + O(-) (17) 
2 + 4r 2 a or 



where r = 4p a /3f3f, and the coefficients are 



3 n 9, , (28 + 157r + 102r 2 ) 
co = --In7r+-ln3-31n2 + ^ - 1 + 3K 2 , (18) 



68 + 1115r + 804r 2 + 1692r 3 9(3 + r) 
Cl ~ 32(1 + 2r)3 4(1 + 2r) (19) 



and 



» .0798. (21) 

It is assuring to note that the coefficient of the 0(1/ a) term is ~ O(l) - large contributions 

from Fi_i oop and F 2 -i oop cancel to give a well-behaved series. In the following, we will use 
eqn(|T7|) for the free energy in this region. 



Pure Yang-Mills theory has other saddle point solutions[12] : the "fluxon" configurations, 
which are topological excitations corresponding to center of the gauge group. These are stable 
for p a > @. For the region of large j3f, contributions of these maxima are much smaller 



than the higher loop terms and so we will neglect these configurations. 

Since the term proportional to (3 a in (|l]) is blind to the center of the group, such configu- 
rations are not suppressed on or near the /3f = axis, and will have to be taken into account. 
This is done || by expanding the partition function in a series in (3f, to obtain 

FmiPfrfa) = "6/3/(1 - \f3 f ) + F som ((3 a + + 0(pj) (22) 
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where, up to one loop, 

F SO (3)(Pa) « 3(p a (v 2 +v 4 -2) + ln(/ (a) - h(a)) + | - 

+- (in 2 + 3 ln(l + 2v 2 ) - 2 ln(l + 3v 2 ) + 2K^. (23) 



Here iTi is given by eqn(|l4|) with a substitution of c = in (15). a,v are solutions of 

lh(a)-I 2 {a) 
3 J (a) -Ii(a)' 



a 



-/3 a w(l + 2v 2 ). 



(24) 
(25) 



One can check that coefficient of the term proportional to 1/a in F so ^ is ~ O(l) - the 
loop expansion does not differ substantially from an expansion in 1/a. Figs. 1 and 2 display 
a comparison of the Monte Carlo data p[ || with the predictions for phase transition lines 
obtained by comparing the free energies Fj, Fji, Fjjjin the different regions. 



Pa 




Fig.l : The points (joined by thin dots) are the Monte Carlo data of Q, with open diamonds 
denoting the points where a first order bulk phase transition was ruled out by |J . The full line 
is the curve obtained in [||, thick dotted line is obtained taking upto 10th order in eqn(||). 



In Fig.l the Monte Carlo data are shown along with the curve obtained in || by comparing 



Fi upto 10th order in eqn(g)and Fu,Fiu upto one loop order in eqns (11), (22). Also shown 
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is the curve obtained by taking the strong coupling series upto 16th order. While taking the 
higher order terms in strong coupling expansion changes the location of the endpoint of the 
phase transition line drastically, moving it closer to the Monte Carlo data, the location of the 
transition line itself changes very little for (5 a > 1.8. 
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Fig. 2 : The Monte Carlo data are shown against prediction using 16th order strong coupling 
expansion and eqns (|l7|) , (p2[) (full line). The thick dotted line corresponds to 10th order in 
strong coupling series. 

In Fig. 2 we show the prediction for phase transition curve obtained by comparing the 
improved mean field series of eqn([T7|), the strong coupling series upto 16th order and eqn(f22|). 
Also shown is the curve for taking strong coupling series upto 10th order. The same feature, 
namely, the extreme sensitivity of the endpoint and robustness of the upper part of the curve 
is noted. Also, a comparison of Figs. 1 and 2 reveals that the improvement in the mean field 
series leads to a curve that is considerably closer to the Monte Carlo results. The convergence 
of the mean field series was checked by checking that on taking just the first two terms in 
(|l7|) , the transition line is left unchanged. The endpoint of the transition line, however, is still 
in disagreement with ||. Due to the extreme sensitivity of the endpoint, it seems that very 
high orders in strong coupling series may be needed in order to pinpoint it precisely. However, 
since the strong coupling expansion seems to be well behaved near the phase transition line for 
Pa > 1-5, the existence of the bulk transition at least upto this region seems to be confirmed. 
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The only caveat for this conclusion is the choice of gauge fixing. It is always advisable to 
check for possible gauge dependence of a result obtained by gauge fixing. In this case, it is even 
more so since the axial gauge constrains the deconfinement order parameter to be nonzero in 
region II whereas it is zero (or small) in the strong coupling region I. It would be interesting 
to confirm the bulk phase diagram in Fig. 2 by using another gauge condition to establish the 
bulk transition for (3 a ~ 1.5 beyond doubt. 
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